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Abstract. The structure and stability of strongly charged complex fibers, formed by complexation of a 
single long semi-flexible polyelectrolyte chain and many oppositely charged spherical macroions, are inves- 
tigated numerically at the ground-state level using a chain-sphere cell model. The model takes into account 
chain elasticity as well as electrostatic interactions between charged spheres and chain segments. Using a 
numerical optimization method based on a periodically repeated unit cell, we obtain fiber configurations 
that minimize the total energy. The optimal fiber configurations exhibit a variety of helical structures for 
the arrangement of macroions including zig-zag, solenoidal and beads-on-a-string patterns. These struc- 
tures are determined by a competition between attraction between spheres and the polyelectrolyte chain 
(which favors chain wrapping around the spheres) , chain bending and electrostatic repulsion between chain 
segments (which favor unwrapping of the chain), and the interactions between neighboring sphere-chain 
complexes which can be attractive or repulsive depending on the system parameters such as medium salt 
concentration, macroion charge and chain length per macroion (linker size). At about physiological salt 
concentration, dense zig-zag patterns are found to be energetically most stable when parameters appro- 
priate for the DNA-histone system in a chromatin fiber are adopted. In fact, the predicted fiber diameter 
in this regime is found to be around 30 nanometers, which appears to agree with the thickness observed 
in in vitro experiments on chromatin. We also find a macroion (histone) density of 5-6 per 11 nm which 
agrees with the zig-zag or cross-linker models of chromatin. Since our study deals primarily with a generic 
chain-sphere model, these findings suggest that structures similar to those found for chromatin should 
also be observable for polyelectrolyte- macroion complexes formed in solutions of DNA and synthetic nano- 
colloids of opposite charge. In the ensemble where the mean linear density of spheres on the chain is fixed, 
the present model predicts a phase separation at intermediate salt concentrations into a densely packed 
complex phase and a dilute phase. 

PACS. 87.15.-v Biomolecules: structure and physical properties - 87.16.Sr Chromosomes, histones - 82.70.- 
y Disperse systems; complex fluids 



1 Introduction 

Recent years have witnessed a growing interest in the 
structure and phase behavior of large charged polymer- 
macroion clusters, for instance, a long polyelectrolyte (PE) 
chain is complexed with many oppositely charged spheres 

striking example is realized in biology, where a huge fiber 
is formed by complexation of a long (negatively charged) 
DNA chain with (positively charged) histone proteins, giv- 
ing rise to the prominent chromatin fiber 16, 17,18, 19, 20, 
I2"l"ir2"2"ll2"5] . The fundamental unit of chromatin is known as 
the nucleosome consisting of about 200 base pairs (about 
68 nm) of DNA associated with a cylindrical- wedge-shaped 
histone octamer of diameter about 7 nm and mean height 
5.5 nm. The nucleosome has two main parts: i) the nucle- 
osome core particle (or the core particle) comprising 146 



base pairs of DNA wrapped in nearly a l-and-3/4 left- 
handed helical turn around the histone octamer, and ii) 
the linker DNA, which connects adjacent core particles to 
one another. There is an additional component known as 
linker histone HI (or its variant H5), which binds to the 
DNA and the histone octamer in such a way that it brings 
the entering and exiting strands of DNA together along a 
short distance, forming a stem-like structure. The role of 
linker histone will be discussed further below. 

In vitro experiments have revealed striking salt-induced 
conformational changes for the chromatin fiber (see, e.g., 
|16[|24[|25[|26ll2"Tll28[|29j ) : At very small salt concentration, 
the fiber structure is swollen, displaying an open beads- 
on-a-string pattern in which individual core particles be- 
come highly separated from each other; the resultant fiber 
is known as the 10 nm fiber because of the nearly 10 nm 
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hard-core diameter of the core particle (nearly equal to the tem throughout this study. We should emphasize, how- 
combined diameter of the core histone octamer and twice ever, that the present model is not meant to address the 
that of the DNA wrapped around it). As the salt concen- specific problem of the chromatin fiber (as no attempt is 
tration increases, the fiber becomes more and more folded made in order to incorporate explicitly the specific struc- 
and finally within the physiologically relevant regime (aroundtural details of the actual chromatin fiber such as the 



100 mM NaCl), it exhibits a thick and dense fiber of di- 
ameter about 30 nm, known as the 30 nm fiber. This salt- 
dependent behavior indicates that predominant electro- 
static mechanisms are involved and influence the organi- 
zation of the chromatin fiber. 

The precise arrangement of nucleosome core particles 
in the 30 nm fiber has been under debate for a long time 



linker histone, core histone tails, specific binding effects 
or the precise shape of the histone octamer [I6lll8j ). but 
we rather focus on the more general aspects of complex 
fiber formation in PE-macroion systems. The results are 
thus also relevant for mixtures of DNA and oppositely 
charged nano-colloids or synthetic spheres. 

We show that the transition between the predicted he- 
[T71ITgll2lll!r71l251l2^^ patterns is dictated by the wrapping-dewrapping be- 



43 ,44 , 45 , 46 , 47 , 48 , 49 , 50 , 5 1 , 52 , 53l . A large number of mod- 
els have been proposed for the structure of the 30 nm fiber 
which are more or less consistent with the experimental 
observations. For instance, several studies [27,28,29,30 
appear to support a zig-zag or cross-linker pattern, where 
successive nucleosomes are at opposite sides of the fiber 
and the linker DNA joining them bridges across the he- 
lical core of the fiber. Among other models supported by 
recent experiments are the solenoidal model [2~4"1[3"T] (in 
which successive nucleosomes are juxtaposed as nearest 
neighbors with the linker DNA bent in between them), 
the twisted ribbon model [32] . as well as the more gen- 
eral class of ribbon models [SUEE] (see Refs. Q1IM1MI] 
for a more comprehensive reference list and discussion of 
various chromatin models). 

Our goal in this work is to address structural prop- 
erties of generic chain-sphere complex fibers comprising 
a long semi-flexible charged PE chain and many oppo- 
sitely charged spheres. We adopt a simple chain-sphere 
cell-model approach [541155] that allows for a systematic 
description of the polymer conformation in an infinite 
fiber. The chain conformation and the fiber structure are 
determined by a competition between chain elasticity, elec- 
trostatic interactions (among chain segments and charged 
spheres across a salt solution) as well as geometrical con- 
straints. We employ a numerical optimization method to 
find configurations minimizing the total fiber Hamilto- 
nian with respect to the whole conformation of the com- 
plexed chain and the position of macroions. This scheme 
is known as the ground- state- dominance approximation 
which is valid for strongly-coupled charged complexes (such 
as DNA-histone complexes) that exhibit large PE adsorp- 
tion energy and relatively small thermal fluctuations (541 
155) . The organization of the chain on individual macroions 
and hence throughout the fiber is not externally imposed, 
but rather obtained as a result of the interplay between 
inter- and infra-fiber interactions. We demonstrate that 
the resultant optimal fiber structures exhibit a variety 
of helical structures including, for instance, beads-on-a- 
string, zig-zag, and solenoidal patterns depending on a few 
basic system parameters, namely, the salt concentration, 
the macroion charge valency, and the chain length per 
macroion (or the linker length). Other parameters such as 
the chain persistence length and its charge density as well 
as the sphere diameter will be fixed; here we choose the 
parameter values appropriate for the DNA-histone sys- 



havior of the chain as the system parameters vary. For 
instance, at low salt concentration, the dominant self- 
repulsion of chain segments leads to dewrapping of the 
chain from macroions, generating an expanded fiber. In 
the physiological regime, zig-zag patterns are found to be 
energetically most stable, because the chain is wrapped 
around macroions in about a l-and-3/4 turn for moderate 
sphere charge |54j . which turns out to be about the same 
degree of DNA- wrapping in the nucleosome core particles. 
At elevated salt (beyond lOOmM monovalent salt), the 
highly wrapped state of the chain leads to very compact 
solenoidal structures for the fiber. The predicted struc- 
tures will be summarized by presenting a two-angle struc- 
tural diagram in analogy with the two-angle models in 
Refs. [T71I27] . An important result follows from our model 
when the chain (linker) length per macroion is treated as 
a free parameter. It is found that the binding energy (per 
macroion) takes a non-convex shape at intermediate salt 
concentrations (when plotted as a function of the chain 
length per macroion), with a global minimum at a high 
sphere density, indicating a gas-liquid-type "phase coexis- 
tence" along the fiber: Part of the PE chain forms a dense 
complex fiber with macroions phase-separating from a di- 
lute phase along the fiber consisting only of uncomplexed 
chain. 

In the aforementioned optimization model, the wrap- 
ping structure of the chain around individual macroions 
is not fixed and can vary so as to minimize the effec- 
tive Hamiltonian of the fiber (unconstrained optimiza- 
tion) . This can be thought of as a simple model mimicking 
linker- histone-depleted chromatin [53]. The 30 nm fiber 
obtained within the unconstrained optimization model shows 
a rather small density (number per projected unit length) 
of macroions along the fiber, i.e., about one macroion per 
11 nm, which is comparable with the nucleosome density 
at relatively small salt concentration (where the chromatin 
still has a 30 nm diameter [H]): but not with the density of 
about 6 histones per 11 nm P^l[2^1[2^1[271[^8l [2^] observed 
in the physiological salt regime. 

As noted above, the linker histone HI "glues" together 
the entering and exiting strands of the DNA and also sta- 
bilizes ("seals off") the internal structure of the DNA in 
the nucleosome core particle. It plays a key role in the com- 
paction of the chromatin into a 30 nm fiber [24] . In order 
to mimic this kind of effect in a very simple way within 
the present model, we also employ a variant of the opti- 
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Fig. 1. Cell model for a complex fiber: a) Each unit cell consists of a single PE-macroion complex the precise configuration 
of which is obtained from a numerical minimization method, b) Two adjacent unit cells (which are identical in configuration) 
are linked together such that they share a common tangent vector at the unit cells junction (see text), c) By repeating this 
procedure the whole complex fiber is constructed from a main unit cell. Within the full minimization scheme, all degrees of 
freedom associated with a single unit cell (that is the chain beads locations as well as the boundary rotational angle x) are 
varied in order to find the optimal configuration minimizing the full effective Hamiltonian of the fiber. 



mization scheme, the so-called constrained optimization, 
in which the core complex structure is fixed according to 
the optimal structure of an isolated complex (with a chain 
length of 146 DNA base pairs); then only the relative ori- 
entation angle of adjacent core particle complexes can vary 
in order to find the optimal fiber structure. This leads to 
a number of significant changes. Most notably, the enter- 
ing and exiting strands of the chain tend to form a cross 
pattern at physiological salt (producing a small entry-exit 
angle as compared with the unconstrained model), which 
resembles the stem-like pattern proposed for Hl-stabilizcd 
nucleosome 24,27,28 . As a result, the resultant optimal 
fiber is found to have a more compact zig-zag pattern with 
a histone density of about 5-6 histones per 11 nm in agree- 
ment with Refs. fHtRMMMM- 

The organization of this paper is as follows: In Section 
[3J we describe the theoretical model and the numerical 
method used to study infinite charged complex fibers. In 
Section [3l the salt-induced structural changes of the fiber 
are studied within the unconstrained optimization model. 
The constrained optimization model is analyzed in the fol- 
lowing Sectional We then show in Section [5] how the fiber 
structure is influenced by variation of the linker length 
and discuss different thermodynamic ensembles and the 
possibility of phase separation along the chain. 



2 Theoretical model: Cell-model approach 

The model used here for the study of complex fibers is 
based upon the chain-sphere model [53]!55] , which consists 
of a single semi-flexible PE chain complexed with an oppo- 
sitely charged sphere (macroion). Here we use this model 
to describe an infinite complex fiber made from complex- 
ation of a single long PE chain with many macroions. We 
restrict the discussion to strongly coupled complexes and 
thus employ the ground-state-dominance approximation 
corresponding to strong adsorption of the PE chain on 
each sphere [54ll55ll57j . 

The main idea here is to apply a cell-model approach 
by describing the complex fiber in terms of identical chain- 
sphere unit cells linked together via the same PE strand. 



In analogy with the chromatin fiber, one can thus think of 
each unit cell as consisting of two main parts: i) the core- 
particle complex (macroion sphere complexed with a chain 
segment), and ii) the linker chain connecting adjacent core 
particles as shown in Fig. [TJ 

In our numerical analysis, the complex fiber is con- 
structed by replicating one unit cell (referred to as the 
main unit cell) infinitely many times such that the con- 
nectivity and smoothness of the linker chain is preserved. 
We begin the minimization with a given initial condition, 
typically, the one in which the structure of the core parti- 
cle complex in each unit cell corresponds to the energet- 
ically optimized configuration in an isolated chain-sphere 
complex with the chain length associated with each core 
particle being taken to be equivalent to 146 DNA base 
pairs. This procedure for a single isolated complex has 
been studied thoroughly in our previous works [54,55,56 . 
We then connect an arbitrary length of linker chain and 
the interaction Hamiltonian of the whole fiber is minimized 
using a recursive numerical method to find the optimal 
fiber structure without any further constraints (see be- 
low). This follows simply by focusing on a single unit cell. 
Although this prescribes a model consisting of identical 
unit cells, the optimized fiber structures are by no means 
trivial in the way the unit cells or core particles are ar- 
ranged. They display a rich structural phase diagram of 
helical patterns. Note that these helical patterns may or 
may not exhibit periodic structures (i.e., structures with 
discrete translational symmetry). In fact, the energetically 
optimal structures are merely quasi-periodic, i.e., they do 
not show exact (rational) periodicity but a distribution of 
periodicities with a sharp peak around a mean value of 
periodicity (see Ref. [55] for details). We also note that 
electrostatic interactions between all unit cells are fully 
taken into account, and thus on this level, the only ap- 
proximation of the present cell model is that structural 
symmetry breaking between the unit cells is neglected. 

2.1 Main unit cell 

Each unit cell consists of a piece of PE chain of length 
L c (or an equivalent of Nbp = £ c /(0.34 nm) DNA base 
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pairs each of length 0.34 nm) with linear charge density of 
— r (in units of the elementary charge e) and bare (un- 
screened) mechanical persistence length of £ p . The PE 
chain interacts with a uniformly charged sphere of radius 
R s and charge valency Z in a salt solution of concentration 
C s . Although the present model is quite general, in making 
explicit calculations, we set the parameters appropriate 
for the DNA-histone system by choosing r = 5.88 nm -1 
(maximum dissociation of two elementary charges per base 
pair [58.), l v = 30 nm [54l[59ll60] and R s = 5 nm, while 
L c , Z and C s are varied. Note that the sphere radius of 
R s — 5 nm represents the closest approach of DNA of ra- 
dius 1 nm on a histone octamer of mean radius of about 
4 nm. As in previous studies 54)55], we employ a dis- 
cretization scheme in order to parametrize the chain con- 
formation as required for numerical analysis. The chain is 
discretized using N + 1 discretization points or beads of 
valency q = tL c /(N + 1) located at positions {r°} (with 
i = 0, 1, . . . , N and the super- index denoting the main 
unit cell). This amounts to N rigid chain subunits within 
a unit cell each described by a bond vector u° = r ° — r°_ 1 
of length |u?| = A and two polar and azimuthal angles 9i 
and <j>i with respect to a fixed orthogonal reference [54 ,55 
(N has in general no connection with the number of ac- 
tual monomers) . The location of the central macroion R° 
in the main unit cell is chosen as the origin, i.e., R° = 0. 
We conventionally discretize each natural DNA base pair 
into two discrete subunits; therefore N — 2Ab p and we 
have q = 1 (i.e., one elementary charge per subunit or 
two per base pair) for large N when DNA linear charge 
density is adopted. 



(and in fact the whole fc-th unit cell as a rigid body) 
such that falls onto u^" 1 to fulfill the smoothness of the 
chain at the junction. We then allow the image unit cell fc 
to rotate around this common bond vector by an arbitrary 
(azimuthal) angle Xi see Fig. Q] for an illustration. This 
angle is indeed a rotational degree of freedom for each 
unit cell as a whole in addition to other intra-cell degrees 
of freedom (bead and sphere positions) that is varied in 
order to minimize the fiber Hamiltonian (see Section 12.41 
below). 

The rotation matrix E is obtained as 

E = BX A (3) 
with the matrices A, X and B defined as follows 



A = 



X = 



COS (pi 

— sin 4>\ cos 9i sin ( 
sin (pi sin 9\ cost 



B 



— sin 0i 

— cos (pi cos 9i 
sin 9 1 cos <pi 

' cos x — sin X 
sin% cosx 
1 

— sin <pjy — cos (pM cos 9m sin 9m cos (pM 
cos (pM — sin <pM cos 9n sin tpN sin^jv 

sin 9 m cos 9 m 



(4) 
(5) 
(6) 



(pi are the polar and azimuthal angles of ii® 
m are the polar and azimuthal angles of u N . 



where 6\ and 
and 9 m and 4 

Note that the matrix E is only a function of the angles 
{x,9i,<f>i,6 N ,(p N }. 



2.2 Unit cell replication 

Let us assume that the structure of the core-particle com- 
plex within the main unit cell is given (from a minimiza- 
tion scheme as discussed later) . We construct the complex 
fiber by attaching an image of the main unit cell to it such 
that the linker chain tangentially connects the two unit 
cells with no kinks at the junction (see Fig. [T]); i.e., the 
last chain bond vector, u N , in the main unit cell is chosen 
as the first bond vector, u}, of the first immediate im- 
age cell (labeled by super-index 1) and so on for all other 
image cells in a recursive manner. 

Formally, this replication procedure is established by 
a combined translation and rotation transformation. In 
general, the position of the i-th chain bead and that of 
the sphere center in the fc-th unit cell (labeled by super- 
index fc), i.e., r^" and R fc respectively, are obtained as 

r^E-^-rg-^+r^ (1) 

R* = E.(R*- 1 -r*- 1 ) + r* r -_ 1 1 , (2) 

from the given configurations in the (fc — l)-th unit cell 
(labeled by super- index fc — 1). Here E is the full rotation 
matrix around the common point or junction r5T i = i"q °f 
the two adjacent unit cells, which rotates the bond vector 



2.3 Intra-fiber interactions 

Because the infinite complex fiber is obtained by repli- 
cating identical unit cells, one can focus on the effective 
Hamiltonian of the main unit cell, which can be written 

as 

M 

U = U se if + J2^ (?) 

k=l 

where H se if gives the self-energy of the main unit cell, and 
Hint accounts for the interactions between the main and 
the fc-th unit cells. The following interactions are included 
in these two terms: 1) electrostatic interactions between 
all charged beads along the PE chain with one another 
and with charged spheres within the linear Debye-Hiickel 
theory (which effectively captures salt screening effects) 
as well as the mutual interactions between the charged 
spheres themselves, 2) a short-ranged soft-core excluded- 
volume repulsion between chain beads and spheres only, 
which prevents the chain from penetrating the spheres, 
and 3) the mechanical bending elasticity of the PE chain. 
Note that the electrostatic contribution to the bending 
rigidity is explicitly included via the electrostatic self- 
energy term. In the second term in Eq. ([7]), we have used 
a cut-off M to account for the interaction between M con- 
secutive unit cells. We typically choose M = 5. As explic- 
itly checked, the results are not influenced by increasing M 
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£ N ( 

Haeli = ^ < 1 - cos(9i — 9i--i) + sm9i sin6>;_i ( 1 - cos(0i — 

i=2 



(8) 



+ q 



i=0 j=i+l 



Zq£ B 
1 + KR* 



E 



- K (\r°\-Rs) 



Ae 



-(|r°|-i? s )/« 



rtint — 1 



B Z^Z^ Irp-r*! 1 + Kj R s L~, 



i=0 j=2 



N 

+£ 

i=2 



„-K(|r fc |-H B ) 

6 ,, - -Ae-d^-^ 



+ 



<|r°-R fc |-fl 3 ) 

|r? -R fc 



(1 + kR s ) 2 |R fc | 



4 e -(|r°-R fc |-ils)/« 



(9) 



further due to the highly screened nature of electrostatic 
interactions at high salt and/or highly expanded size of 
unit cells at low salt. 

The self-energy of the main unit cell is written in dis- 
cretized form as 

See equation (0) above 

(in units of fee? 1 ), where the first term is the mechani- 
cal bending contribution, the second term is the electro- 
static repulsion between chain segments, and the third 
term represents chain-sphere attraction as well as the soft- 
core repulsion (which is specified by a repulsion range, 
a, and strength A, which are fixed at typical values of 
A = 0.014 nm" 1 and a = 0.02 nm 56 ). Recall that 
k = (8tt£bC s ) 1 ^ 2 is the inverse Debye screening length 
for monovalent salt with 1-q — e 2 / {Anee^k-QT) being the 
Bjerrum length. 

The interaction between the main unit cell and the 
fe-th image unit cell reads (in units of k^T) 

See equation above 

noting that the first two beads i — 0, 1 for the fe-th image 
unit cell are counted as the last beads of the (k — l)-th 
cell as explained before. The first term in Eq. (|9]) corre- 
sponds to repulsion between chain segments in the main 
unit cell and those in the fe-th image cell, the second 
term involves electrostatic attraction and soft-repulsion 
between the chain segment in the main unit cell and the 
sphere in the fe-th image cell and vice versa. The last term 
represents electrostatic repulsion between the spheres in 
the two unit cells. 



2.4 Numerical method: Unconstrained (full) 
minimization model 

For strongly coupled PE-macroion complexes (i.e., with 
large polymer adsorption energy), the so-called ground- 
state-dominance approximation becomes accurate [54 . Here, 
one seeks the optimal or ground-state configuration of the 
fiber minimizing the effective Hamiltonian ((7} with respect 
to all degrees of freedom that are necessary to specify the 



overall chain conformation. This amounts to 2N + 1 inde- 
pendent degrees of freedom (per unit cell) {9i, <fti=£p, rp, x}, 
where 9i and 4>i^p are the chain's bond vector angles, rp is 
the z component of the position of the middle bead of the 
chain P = N/2 (measured from the sphere center) and 
X the unit cell rotation angle. Note that 4>p (azimuthal 
angle of the subunit P or bond vector up) is fixed in or- 
der to remove Goldestone (zero) modes associated with a 
trivial rotational symmetry around the z axis connecting 
the sphere center to the middle bead P. 

Due to a large number of free variables, we tackle the 
problem numerically using the quasi- Newton optimization 
algorithm |61] , which locates the minimum of the Hamil- 
tonian using its first-order derivatives. At a minimum, the 
first-order variation of the Hamiltonian vanishes, i.e. 

Sn = >f -rrS9i+— 5fa P )+^X+^Srp = 0. 

^ \d6i d(j)^p ) dx drp 

( 10 ) 

We use a combination of stochastic and parameter 
quenching methods to make sure that the obtained struc- 
tures are not meta-stable (local minima) . The meta-stable 
structures have been discussed elsewhere for chain-sphere 
complexes [541I56] . 

3 Salt-induced structural changes of the 
fiber - unconstrained optimization 

In this section, we investigate the dependence of structural 
properties of the complex fiber on the salt concentration 
for fixed macroion charge valency Z = 15 and fixed to- 
tal length of chain per unit cell L c = 68 nm (equiva- 
lently 7V bp = 200 base pairs of DNA) [251128]. The salt 
concentration C s is increased from zero up to the physi- 
ologically relevant range of 100 mM monovalent salt and 
beyond (with the inverse screening length spanning the 
range k = — 1.1 nm -1 ). Here we focus on the uncon- 
strained (full) minimization scheme; therefore, the chain 
length wrapped around the macroion is not fixed and can 
change accordingly. The results display a variety of helical 
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Fig. 2. The optimal spatial configuration of the complex fiber as calculated from the unconstrained minimization of the 
Hamiltonian Eq. (J7J) for macroion (sphere) charge valency Z — 15, total chain length per unit cell (per sphere) of L c — 68 nm 
(equivalent to 200 DNA base pairs) and for various Debye inverse screening lengths as indicated on the graph. 



fibers as visualized in Fig. [2J which reflect drastic struc- 
tural variations by changing the salt concentration. 

The ground-state configurations shown in the Figure 
can be grouped into a few generic classes, such as beads-on- 
a-string (Figs, [^i, b, g), zig-zag or cross-linker (Figs. [^Ji, 
h-1) and loose solenoidal structures (Figs.[2j:, e, f). These 
types of patterns have been identified more specifically 
within recent geometric models for the chromatin fiber 18 , 
|2"7II35| . In general, the zig-zag pattern can be thought of 
as a configuration in which successive core particles are lo- 
cated on opposite sides of the fiber and the straight linker 
chain joining them bridges across the helical core of the 
fiber. This property is absent in the loose solenoidal pat- 
tern where successive nucleosomes are adjacent (as nearest 
neighbors) . Note that this classification is only qualitative 
in the present context and will be used for convenience 
when referring to the structures obtained from our model. 
A more quantitative analysis of the salt-dependent behav- 
ior will be discussed in the following Sections. 

As seen, the degree of chain complexation increases 
with increasing salt concentration. The highly wrapped 
states occur at intermediate salt concentration and result 
from the growing dominance of electrostatic chain-sphere 
attraction against the chain (electrostatic and bending) 
self-energy. Indeed, as we show later, macroions are highly 
overcharged by the PE adsorption much in the same way 
as a single isolated complex is overcharged at intermediate 
to large salt concentration 54,55,62,63 . 



At vanishing salt concentration (Fig. [5^), the PE chain 
takes a straight-line shape due to the dominant electro- 
static self-energy of the PE chain as charged units inter- 
act with long-ranged (unscreened) Coulomb interaction. 
The distribution of spheres around the chain is periodic 
but their centers lie on different angular locations around 
the chain. For inverse screening length in the range of 
k = 0.3 to 0.6 nm" 1 (Figs.[5Ji-g), the fiber structure shows 
an anomalous behavior, i.e., the diameter of the fiber de- 
creases and the projected distance (along the fiber axis) 
between neighboring sphere centers increases, which is due 
to the fact that the chain strand completes its first turn 
around each sphere (see also Fig. [3]). Upon further com- 
plexation of the chain (Figs.[Hi-l), the fiber thickness tends 
to increase again and the spheres become more densely 
packed. Note that in the range of Debye inverse screen- 
ing lengths k — — 0.5 nm -1 (<25 mM monovalent salt), 
the chain wrapping is less than a complete turn and these 
configurations appear locally as beads-on-a-string struc- 
tures. However, the chain wrapping may be increased by 
taking a higher sphere charge, leading to different struc- 
tures even for small shown elsewhere |5BJ. The role 
of sphere charge will be discussed later. 

3.1 Overcharging of the core particle 

It is useful to consider the local structure of the core par- 
ticles within the fiber in more detail. A closer view of the 
unit cell of the fiber is shown in Fig. [3] for different salt 
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Fig. 3. The local shape of the core particles within the optimal fiber structures in Fig. [2] for sphere charge valency Z = 15 and 
various Debye inverse screening lengths as indicated on the graph. The total length of the PE chain per unit cell is fixed as 
L c — 68 nm (JVb p = 200). The increase in the number of chain turns around the macroion is clearly demonstrated. 
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Fig. 4. Number of chain beads (circles) that are attached to 
the sphere (i.e., within the distance 1.02J? S from the sphere 
center where R s = 5 nm is the sphere radius) as a function of 
the inverse screening length for the structures shown in Figs. 
[2] and [3] As noted before each bead is monovalent (q — 1), 
thus this number equals the PE charge adsorbed on the sphere. 
Diamonds show the net charge of the complex. At very low salt 
concentration, the sphere charge is already nearly compensated 
by the adsorbed PE segment and for larger salt, the sphere is 
strongly overcharged. 



almost linearly with the salt concentration and reaches 
a plateau-like region at high salt in agreement with the 
behavior observed for a single isolated complex [56]. The 
saturation at high salt is due to the finite and fixed length 
of the PE chain per sphere that is 200 DNA base pairs 
in the present case. The adsorbed PE charge amounts to 
that of a DNA segment of 146 base pairs at about the 
physiological salt concentration (k ~ 1 nm -1 ) 54J . 

Note that for the chosen parameters here (with sphere 
valency of Z = 15), the core particle is always overcharged, 
that is U — Z > 0. The overcharging degree, (U — Z)/Z, 
becomes quite large at intermediate salt, such that the net 
complex charge becomes more than twenty times larger in 
magnitude than the bare sphere charge. This striking fea- 
ture is known to occur due to strong lateral chain segment 
correlations at large salt concentration, which leads to a 
highly ordered adsorbed layer on the macroion surface [SH 
55,62,63, 64 . In this situation, the spacing between chain 
strands on the sphere is of the order of the Debye screening 
length, k,- 1 [53]. 

The unit cell configurations in Fig. [3] also display how 
the relative entry-exit angle of the chain in a core parti- 
cle is influenced by the changes in the salt concentration. 
The changes in the relative entry-exit angle in turn affects 
the overall arrangement of the core particles and thus the 
global structure of the fiber as seen in Fig. [5] 



concentrations. We introduce two quantities that measure 
the degree of PE chain wrapping around each macroion 
in the fiber, namely, the adsorbed PE charge, U (that in- 
cludes only those chain beads that are within distances 
less than or equal to 1.02i? s from the sphere center), and 
the net charge of the core particle defined as U — Z. Both 
quantities are shown in Fig. [4] as a function of k. 

The unit cell configurations clearly show that the PE 
chain becomes gradually more wrapped around the sphere 
upon increasing the salt concentration. That amounts to 
more than two complete turns at high salt (k — 1.5 nm -1 ). 
As seen from Fig. [H the adsorbed PE charge increases 



3.2 Geometry of the fiber 

Besides the overall shape, one is interested in geometrical 
characteristics of the fiber and the way they vary with the 
salt concentration. Here we consider the following param- 
eters (see Fig. [5]): 

- The fiber diameter, Dp, which refers to the diameter 
of the smallest outer cylinder coaxially enclosing the 
fiber. 

- The projected center-to-center distance, du, of consec- 
utive spheres along the central fiber axis, which gives 
a measure of projected density of the fiber, n\\; conven- 
tionally, we represent the fiber density as the number 
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Fig. 5. Schematic representation of the fiber diameter, Dp, and the projected center-to-center distance, dn, of consecutive 
spheres along the fiber axis. The central fiber axis is shown by a horizontal solid line. 
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Fig. 6. a) The projected distance of neighboring sphere centers along the fiber axis, d\\, and b) the fiber diameter, Df, as a 
function of the inverse screening length. The macroion charge and chain length per unit cell are fixed as Z = 15 and L c — 68 nm. 



of spheres per 11 nm as typically expressed in experi- 
mental literature for chromatin I28j : i.e. 



where du is measured in units of nm. 
- The entry-exit angle, ip, and the dihedral angle, £, 
which can be defined (Section @| in analogy with the 
chromatin models (see Fig. [7]). 

As seen in Fig.|BJi, the projected distance takes a max- 
imum value of about du =68 nm at zero salt concentra- 
tion and then falls off in a non-monotonic fashion to quite 
small values of about sphere radius du — R s = 5 nm, when 
the inverse screening length (salt concentration) exceeds 
k = 1 nm -1 (100 mM NaCl) representing a dense fiber of 
a density about nu = 2. The projected distance displays a 
local minimum and a local maximum at inverse screening 
lengths of about n = 0.3 nm -1 and k = 0.6 nm" 1 , re- 
spectively. These extrema roughly correspond to extrema 
in the fiber diameter, D-p, which shows an almost oppo- 
site trend as compared with the projected distance when 
the salt concentration is varied (Fig.[B}3). At zero salt, the 
outer fiber diameter is roughly twice the sphere diame- 
ter, i.e., Dp — 20 nm, due to a staggered configuration 
of the spheres. It changes rapidly at small salt concen- 



tration with k, reaching a maximum and then a sharp 
minimum at an intermediate salt concentration of about 
k = 0.6 nm , where the projected distance exhibits a lo- 
cal maximum. The minimum fiber diameter at this point 
is about 10 nm (equal to the sphere diameter); hence at 
k = 0.6 nm -1 , the complex fiber shows its largest pro- 
jected distance to diameter ratio du/Dp = 4, reflecting an 
almost one-dimensional beads-on-a-string structure (see 
Pig. IS). 

Interestingly, over the range of intermediate salt con- 
centrations between k = 0.8 and 1.3 ran -1 (correspond- 
ing to 60-160 mM monovalent salt), the fiber diameter 
remains almost unchanged at about Dp = 30 nm (re- 
flected by a plateau- like region in Fig. ®>), although the 
fiber structure changes from the zig-zag to loose solenoidal 
pattern (see also Fig. [8]). The chain wrapping in this re- 
gion amounts to about l-and-3/4 turn, close to that in 
the nucleosome core particles in chromatin at intermedi- 
ate salt concentration about the physiological regime [Pol 
l2"4"ir2"TU2"S] . In fact, the experiments reveal that the chro- 
matin fiber exhibits a dense structure in this regime with 
a diameter of about 30 nm, the so-called "30 nm fiber". 
The present results, which predict a stable 30 nm com- 
plex fiber at intermediate salt concentration, thus appear 
to give a trend consistent with those found in chromatin 
experiments ! 16 , 24 ,27, 25] . The present full-minimization 
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Fig. 7. A schematic view of the two-angle representation of macroion positions along the chromatin fiber. Each two successive 
spherical macroions, labeled by n — 1 and n, are connected to each other via a vector a„. The entry-exit angle ip is defined as 
the complementary angle between any two consecutive connecting vectors (shown in the figure for a n +i and a n +2), which is the 
same for all n in a homogeneous helical structure as is the case in the present model. On the other hand, each three successive 
macroions (e.g., those labeled by n — 1, n and n + 1) define a plane with a normal vector b n = a„ x a„+i. The dihedral angle 
^ is defined as the angle between two such normal vectors (shown in the figure for b n and b„+i). 



model however is very different from the Hl-stabilized 
chromatin because (aside from other realistic factors in 
chromatin that are missing in our model) the structure of 
the core particle (and the conformation of the entering and 
exiting chain strands) can change subject to interactions 
with other units. Thus the local fiber structure is different 
from that of chromatin. We also obtain a lower density of 
spheres (tin < 1) as compared with the chromatin fiber. 
The present results however suggest that such 30 nm fibers 
may be feasible in mixtures containing DNA and synthetic 
nano-colloids at intermediate salt concentrations. 

For large salt concentrations (beyond re — 1.3 nm -1 ), 
the fiber diameter rapidly increases but the projected dis- 
tance between spheres remains constant. At re = 1.5 nm , 
the fiber exhibits a small ratio of dn/D-p — 0.077, corre- 
sponding to a tight solenoidal shape (shown on Fig. [8]). 

3.3 Two-angle structural diagram 

The arrangement of macroions along the fiber within the 
present model may also be described using a two-angle 
description in the spirit of the two-angle model presented 
in Ref. [37] for the chromatin fiber. For the sake of analogy, 
we define an entry-exit angle, ip, and a dihedral angle, £, 
as follows (see Fig. [7J. 

Let us assume that a„ = R™ — R"" 1 is the connecting 
vector between two successive macroion centers, where R n 
is the position of the rt-th sphere. The entry-exit angle, ip, 
is defined from the angle between two such consecutive 
connecting vectors, a„ and a„ +1 , i.e. 

, -1 ( a " ' a "+l \ / 1 r,N 

ip = 7T — cos - — n . (12) 

\|a„||a„ + i|/ 

The dihedral angle, £, the angle by which four consec- 
utive spheres are out-plane, is determined by calculating 
the angle between the normal vectors of two planes, one 
identified by the spheres R" _1 , R™, R™ +1 and the other 
identified by R™, R n+1 , R"+ 2 . Explicitly, one has 

e = cos {K¥^\) (13) 



with h n — a„ x a„ + i. In a homogeneous helical struc- 
ture, which is the case in the present model, ip and £ are 
independent of n and thus the same angles are obtained 
between any two consecutive vectors. 

By calculating these two angles for various salt concen- 
trations, one can sketch a two-parameter £ — 1/> structural 
diagram. The result is shown in Fig. [5] for fixed sphere 
charge Z = 15 and fixed chain length (per unit cell) of 
A^,p = 200. Open circles give the two-angle coordinates 
associated with the numerically calculated optimal con- 
figuration at given inverse screening lengths (indicated by 
numbers on the graph). 

It follows from the diagram that the angle ip deter- 
mines the compactness of the fiber, while £ roughly deter- 
mines its structural class, e.g., being zig-zag or solenoidal. 
At small salt concentrations, these angles vary rapidly 
with the salt concentration, which, as mentioned above, is 
a consequence of the chain wrapping process. At high salt, 
most of the free length of the chain available in a unit cell 
is already wrapped and thus the angles vary weakly. The 
parameter which varies more dramatically with re is ip. As 
seen, ip never becomes smaller than 7r/4, due to excluded- 
volume interactions |18) . while £ spans the whole interval 

< £ < 7T. 

The beads-on-a-string patterns are found at the bound- 
aries of the diagram with entry-exit angles roughly in 
the range 3ir/4 < ip < it. While the compact solenoidal 
structures appear in the lower right corner, that is with 
small dihedral angle < £ < 7r/4 but large entry-exit 
angle 37r/4 < ip < tt. One may thus distinguish com- 
pact solenoidal structures on the lower-right corner of the 
graph. The 30 nm zig-zag structures obtained at interme- 
diate salt regime typically appear on the left upper corner 
of the graph, that is for 7r/4 < ip < n/2 and typically for 
large £ > ir/2. 



4 Constrained optimization model 

As noted before, several features of the actual chromatin 
structure are not accounted for within the present model. 
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Fig. 8. The two-angle diagram for fiber structures from unconstrained optimization at fixed macroion charge valency Z = 15, 
fixed chain length per unit cell L c — 68 nm (Nh P = 200), and for various k. The value of k associated with each point is indicated 
on the graph (in units of nm^ 1 ). The zig-zag structures (upper left corner), beads-on-a-string structures (upper right region), 
and compact solenoids (lower right region) are indicated schematically by closed ellipses and for each case a characteristic 
structure is shown. 




Fig. 9. In the constrained optimization model, the conformation of the chain in each unit cell is fixed and only the rotational 
angle \ can vary. The core particle complex (a) consists of a chain segment of length equivalent to 146 DNA base pairs, 
whose configuration is obtained from numerical minimization of a single isolated complex for a given sphere charge and salt 
concentration. The core particles are linked together with straight linker chain of arbitrary length (b) (in the shown configuration, 
the total chain length per unit cell is 7Vb p = 200). The replication procedure (c and d) is the same as explained in Section ^, 21 



One factor is the so-called linker histone HI, a cationic 
protein which binds near the entry-exit region of the DNA 
in the nucleosome. It brings together the two strands of the 
DNA that enter and exit the central core particle forming 
a stem-like pattern [%4"ll28j . In in vitro experiments, HI hi- 
stone can be removed by exposing the system to high salt 
concentrations (about 0.6 M), and washing the solution 
in order to remove dissociated linker histone. In such Hl- 
depleted cases, the chromatin fiber appears to be rather 
open and loose (i.e., with a low number of histones per 
unit length along the fiber) [24] . While in the experiments 
with Hl-stabilized linker histone, the core particle config- 



uration is fixed and no unwrapping of the DNA on the 
histone octamer is possible [35]. As a result, the number 
of base pairs complexed with the histone octamer is fixed. 

Motivated by this observation, we introduce a vari- 
ant of the preceding model, in which the structure of the 
core chain-sphere complex is constrained to be fixed. At 
any given value of the macroion charge Z and the inverse 
screening length k, we fix the core structure within each 
unit cell according to the optimal (ground-state) configu- 
ration of an isolated complex (with a total chain segment 
of length about 146 DNA base pairs). This structure is 
obtained from a similar minimization study as discussed 
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Fig. 10. Optimal fiber configurations obtained from the constrained optimization model for different sphere charge Z as 
indicated on the graph but for fixed inverse screening length k = 1 nm -1 and PE chain length per unit cell Nbp = 200. The 
inverse projected distance, d^ 1 , the binding energy, AE (in units of ksT), and the fiber diameter Df are given in the graph for 
each fiber. The corresponding local structure of the fiber unit cell is shown on the right. 



in Rcfs. 54,55 . We then connect straight linker chains of 
arbitrary length to the strands of the chain entering and 
exiting the core particle. This constitutes a single unit 
cell for the fiber which is hence constructed by replicating 
this unit cell using the transformation method explained 
in Section I^F see Fig. [SI Note, however, that in this case, 
there is only one degree of freedom, that is the angle %, 
which can be varied in order to find the minimum of the 
effective Hamiltonian of the fiber (see Sections 12.31 and 
12. 4p . Note also that the entry-exit angle is dictated by 
the configuration of a single complex in isolation, which 
thus depends only on the sphere charge, Z, and the salt 
concentration via K, but not on the linker length. 



4.1 Optimal constrained configurations for various 
sphere charge 

In Fig. [TU1 we show several structures obtained from the 
constrained minimization for the inverse screening length 
K = 1 nm -1 and for large to small sphere charge Z (at 
fixed chain length iVb P = 200). It is seen that the compact- 
ness of the fiber is less sensitive to changes in the sphere 



charge. As a result, the entry-exit angle of the PE strand 
changes weakly as compared with the unconstrained min- 
imization case. The fiber becomes denser for large sphere 
charge: For Z = 15. 



we have d, | 1 



= 0.35 nm corre- 



sponding to n\\ — 3.85 spheres per 11 nm length along 
the fiber axis, while for Z = 100, we have nn = 5.5 
(d n 1 = 0.5 nm" 1 ), very close to the experimental chro- 
matin value of about 5-6 histones per 11 nm as reported 
inRefs. [Ml[26l[2a[28j[29]. 

In order to understand the underlying mechanism be- 
hind the compaction of the fiber with increasing the sphere 
charge, we consider the binding energy of the fiber, AE, 
i.e., the energy of the optimal configuration measured with 
respect to the reference state of a free straight PE chain 
in the absence of attached spheres. While for Z = 15, the 
binding energy is positive (which reflects the fact that we 
are performing a constrained optimization which does not 
necessarily find the global energy minimum), it becomes 
highly negative for larger Z as indicated in Fig. [TU1 i.e., 
AE = -580.3 (in units of k B T) for Z = 100. This in fact 
reflects effective attractive interactions between core par- 
ticle complexes within the fiber, which result from correla- 
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Fig. 11. Same as Fig. [10] but for fixed sphere charge Z = 100 and 40 and different inverse screening lengths as indicated on 
the graph. 



tions between positive and negative patches on neighbor- 
ing core complexes as may be seen in the spatial configu- 
rations. Such inter-complex correlations have been studied 
before [55] , 

Note also that the diameter of the fiber is quite stable 
at about Dp = 30 nm, that is about the diameter of the 
Hl-intact 30 nm chromatin fiber in the physiological salt 
regime (i.e., k ~ 1 nm -1 ). In fact, the striking point here 
is that within this simple model, the entering and exiting 
strands of the PE chain form a cross pattern as explicitly 
shown by focusing on a single unit cell in Fig. [TU] (shown 
on the right). This situation-that resembles the experi- 
mentally observed stem-like pattern for the DNA in the 
nucleosome-is only obtained within the constrained mini- 
mization and leads to a small entry-exit angle, and hence 
a rather compact fiber in qualitative agreement with the 
Woodcock model [27] . 



4.2 Optimal constrained configurations for various salt 
concentration 

In Fig. [TTJ we show configurations at fixed Z and different 
salt concentrations. As seen for Z = 100, the structure of 
the fiber becomes expanded for small direct con- 

sequence of the increased electrostatic repulsion between 
the spheres. The density of the fiber reduces to about one 
sphere per 11 nm and the fiber diameter increases to about 
40 nm as indicated on the graph. 

As discussed elsewhere [5"4"ll5"5] , the sphere charge should 
be taken as an effective (or fitting) parameter within the 
present model, because the true sphere charge is regu- 
lated under experimental condition depending on the pH, 
or it may shift from the bare value due to charge renor- 
malization effects in a salt solution [58]. In fact, within the 
present model, the most relevant fiber configurations (that 
might resemble the chromatin fiber best) are obtained by 
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Fig. 12. a) and b) Density of spheres (number of spheres per llnm along the fiber axis) nil, and c) and d) the fiber diameter, 
Df, as a function of the inverse screening length as obtained from the constrained optimization model. The macroion charge is 
Z = 40 for a) and c) and Z = 100 for b) and d). The chain length per unit cell is fixed at 200 DNA base pairs (L c = 68 nm). 



changing Z according to the salt concentration. This may 
be thought of as a fitting procedure, since as explained 
above, the effective sphere charge is expected to depend 
on n. To make this point clear, note that for instance the 
fiber configurations for Z = 100 and low k do not show 
zig-zag patterns [2HI25U28U27U29] . But zig-zag patterns at 
small K may be recovered by taking a smaller value for the 
sphere charge Z (compare k = 0.5 nm -1 for Z = 40 and 
100 in the Figure). 

To summarize these results, we plot the the fiber di- 
ameter, Z?p, and the fiber sphere density, nn, in Fig. [12] 
for two different values of sphere charge Z = 40 and 100. 
Qualitatively, it appears that the sphere density increases 
from about one at low salt up to 5-6 spheres in the physio- 
logical regime of k ~ 1.0 nm -1 (lOOmM monovalent salt) 
[18 , 24 , 2 611271138125] . Examining the fiber diameter in Figs. 
[T2T : and d, we find that the fiber diameter decreases down 
to almost exactly 30 nm at the physiological regime. It 
is noteworthy that the arrangement of the spheres in the 
resultant 30 nm fiber shows a clear double helix pattern 
for Z — 40 and k = 1 nm, -1 . 



5 Variation of the chain length - 
unconstrained optimization 

In this Section, we return to the unconstrained optimiza- 
tion model and address the effects due to changes in the 
length of the PE chain per unit cell, iVb p , on the fiber 
structure. We change Nb p in the range of 50 to 240 DNA 
base pairs and focus on the inverse Debye screening length 



k = 0.8 nm 1 and the sphere charge valency Z = 15. At 
this salt concentration, the fiber exhibits a 30 nm zig-zag 
pattern (Fig.© [53]. 

5.1 Two-angle structural diagram 

In Fig. [T3l we have sketched a two-angle diagram in terms 
of the entry-exit angle, ip, and the dihedral angle, £, sum- 
marizing the structural changes of the fiber with the chain 
length per unit cell, A^p (indicated by numbers on the 
graph). As seen, the optimal configurations are mostly 
gathered in the boundaries of the graph (unlike the be- 
havior found for varying salt concentration in Fig. [5]). 

The global features may be summarized as follows. For 
relatively large chain length per unit cell, the overall struc- 
ture changes weakly with the length and exhibits, for the 
parameters chosen here, the zig-zag pattern (upper left 
corner with large dihedral angles typically £ > ir/2 and 
small to intermediate entry-exit angles 7r/4 < -0 < 57r/8). 
This is because the structure of the core particle is nearly 
independent of the chain length in this regime. Note that 
for iVbp > 200, almost 120 base pairs of the PE chain are 
adsorbed on each sphere, and the rest serves as the linker 
chain. As the chain length decreases down to the length 
adsorbed on the sphere, i.e., Ab p = 120, the diameter of 
the zig-zag fiber reduces to that of the spheres display- 
ing a compact 10 nm fiber with a beads-on-a-string struc- 
ture (upper right corner with ip > 3ir/4 and £ > 37r/4). 
For smaller chain lengths, there is practically no linker 
chain present and the core particles will thus be located 
very close to each other. There is a narrow range of chain 
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Fig. 13. Entry-exit and dihedral angles diagram for various chain lengths per macroion from unconstrained minimization. The 
corresponding values of iVbp are indicated on the graph. The zig-zag region (top left), beads-on-a-string region (top right) and 
compact solenoidal region (bottom right) are shown by closed ellipses. Here re = 0.8 nm -1 and Z — 15 are fixed. 



lengths 50 < Ab P < 100, where the fiber adopts a compact 
solenoidal form, with < £ < ir/2 and ip being distributed 
around 37r/4. Note that the transition between zig-zag and 
beads-on-a-string structures involves an increasing entry- 
exit angle ip which maximizes the distance between adja- 
cent spheres (minimizing their repulsion), while £ changes 
only weakly. In contrast, the transition between beads-on- 
a-string and compact solenoidal structures involves larger 
variations in £. The isolated zig-zag pattern with chain 
length Nbp = 180 corresponds to a local minimum in the 
binding energy of the fiber as will be discussed below. 

5.2 Energetic evidence for intra-fiber phase separation 

In the preceding Section, we discussed configurational vari- 
ations associated with changes of the length of the PE seg- 
ment per unit cell. Now we turn to thermodynamics and 
assume that we have a solution of an infinitely long poly- 
electrolyte and oppositely charged spherical macroions. In 
this case, an interesting question is one of the equilibrium 
state of the system. In other words, what would be the 
resultant linear number density of macroions attached to 
the polyelectrolyte chain and the fiber structure? Clearly, 
a complete answer to this question requires a full thermo- 
dynamic study by considering the free energy in an appro- 
priate ensemble. But within the ground-state-dominance 
approximation, one can address this question on the en- 
ergetic level, i.e., neglecting entropic effects. 



In order to investigate the stability of the system, we 
calculate the fiber binding energy, AE = H — Ho, where 
the first term is the energy of the complex fiber, defined in 
Eq. ([7|) , and the second term is the energy of a free straight 
PE chain in the absence of attached spheres, which is used 
as reference. Note that AE, H and Ho are energies per 
unit cell and in units of fceT. One can then envision two 
different ensembles: An ensemble consisting of a fixed num- 
ber of macroions adsorbed on a fixed large length of the 
PE chain (ensemble A), and an ensemble consisting of 
macroions adsorbed from a bulk reservoir (with a fixed 
chemical potential) on a fixed large length of the PE chain 
(ensemble B). 

Ensemble A is realized in a situation where the total 
number of macroions is finite and thus macroions are de- 
pleted from the solution, and where the total length of 
DNA available for adsorption is fixed such that the mean 
linear concentration of adsorbed macroions is finite and 
can be treated as an input parameter. To study thermo- 
dynamic equilibrium, in this case one has to consider the 
binding energy per unit cell or sphere as a function of the 
chain length per unit cell, A^p- This is shown in Fig. [T5b 
for fixed re = 0.8 nm -1 and Z = 15. As seen, AE appears 
to have two minima at intermediate values of A^p, i.e., 
one at about A^p = 100, where the fiber exhibits a com- 
pact solenoidal structure, and the other minimum at about 
A^bp = 120, where a compact beads-on-a-string structure 
is obtained as discussed before. Also AE shows an over- 
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Fig. 14. Schematic view of the ensemble A with fixed number of macroions adsorbed on a fixed length of the PE chain, 
exhibiting phase-coexistence between a dense and a dilute phase. 




Fig. 15. a) Binding energy of the complex fiber (per unit cell and in units of ksT) as a function of the chain length per unit 
cell iVbp and b) the binding energy per chain unit length as a function of the concentration of adsorbed macroions on the chain, 
1/Nbp, for different values of k as indicated in the graph. 



all non-convex behavior between 7Vb p = 100 and the limit 
of an infinitely dilute fiber (with no macroion adsorbed) , 
i.e., iVbp — > oo. Thermodynamic stability follows from the 
standard common-tangent construction (indicated by the 
broken line in Fig. 115b ). which indicates a gas- liquid- type 
"phase coexistence" along the fiber between a liquid phase 
of high sphere concentration and a gas phase of vanishing 
sphere concentration (see Fig. Il"4"l) . In between, the energy 
has a non-monotonic behavior with values higher than the 
common-tangent line, which thus represents a region of 
meta-stable structures. Adding the one-dimensional trans- 
lational entropy of adsorbed macroions would lead to a 
dilute gas phase with a finite sphere density. Note that 
we do not consider the possibility of a higher density of 
spheres on the chain (i.e. a smaller number N\, p of base 
pairs per unit cell) than corresponding to the minimum in 
energy in Fig. 115a . 

In the case of ensemble B, one deals with a complex 
fiber in equilibrium with a reservoir of macroions, where 
the density of macroions on the chain is controlled by the 
chemical potential of macroions in bulk solution. To study 
thermodynamic equilibrium in this case, one has to con- 
sider the binding energy per unit length of the PE chain 
AE/ A^bp = 'H/A r bp— Wo/Nbp as a function of the adsorbed 
macroion density along the chain, l/A^p, which is shown 



in Fig. H"5b for fixed k = 0.8 nm _1 and Z = 15. We also 
show data for k — 0.4, 1 and 2 nm _1 for large A^p- Note 
that the data corresponding to k = 0.8nm _1 lie below the 
other data, i.e., for this range of A^p, the energy of optimal 
fiber is minimal at n — 0.8nm _1 as compared to other val- 
ues of k [S5]. For the data set for n = 0.8 nm" 1 , AE / 'Ab p 
exhibits two minima at intermediate values of macroion 
density, i.e., one at about 1/ATbp = 0.01 (Ab P = 100), and 
the other minimum at about l/A^p = 0.008 (A^p = 120) 
Note that the location of the minima of AE/N-^p are in 
principle different from the minima of AE, but the relative 
shift is rather negligible. Phase coexistence between two 
fiber phases of different densities in this ensemble would 
correspond in fact to three phase coexistence, i.e. match- 
ing chemical potentials in the two fiber phases and in the 
bulk macroion solution, and will thus be rather unlikely 
since it only occurs when the bulk chemical potential ex- 
actly matches the chemical potentials of the dense and di- 
lute phases simultaneously. In most cases, the equilibrium 
will thus consist of a single phase with a macroion density 
corresponding to the minimum of the energy density in 
Fig. [T5b . If the number of spheres in solution is less than 
needed to homogeneously populate the chain, one crosses 
over to ensemble A with a fixed sphere density which is 
less than the optimal state of ensemble B. 
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6 Conclusion and discussion 

In this paper, we present a systematic numerical approach 
for investigating the structural properties and energetic 
stability of complex hbers formed by complexation of a 
long, semi-flexible polyelectrolyte chain with many oppo- 
sitely charged spheres. The complex fiber is described us- 
ing a chain-sphere cell model, in which the detailed struc- 
ture of the PE chain (locally wrapped around individual 
spheres) as well as interactions between various compo- 
nents (chain segments and spheres) are taken into account. 
These interactions include salt-screened electrostatic in- 
teractions between all charged units (using linear Debye- 
Hiickel potential) , the bending elasticity of the PE chain as 
well as the excluded- volume interactions between the PE 
chain and spheres. Here we choose parameters consistent 
with the DNA-histone system (e.g., using the DNA charge 
and persistence length and a hard-core sphere radius of 
5 nm), but vary the salt concentration, sphere charge, and 
chain length per sphere as control parameters. 

We focus on the optimal structure of the fiber, i.e., 
the structure which minimizes the effective Hamiltonian 
of the system for a given set of parameters. This amounts 
to a ground-state analysis, which is expected to be valid 
for strongly coupled complexes (large polymer adsorption 
energy and relatively small thermal fluctuations) |54[|55] . 
Two different schemes are used here for the optimiza- 
tion procedure: i) unconstrained optimization, and ii) con- 
strained optimization. In the former case, we treat the po- 
sitions of all chain beads within a fiber unit cell as degrees 
of freedom in addition to an independent degree of free- 
dom that defines rotational orientation of two adjacent 
unit cells (core particles) with respect to each other. This 
case therefore requires a many- variable optimization anal- 
ysis. In the other case, the conformation of the PE chain 
on each macroion sphere is fixed according to the optimal 
configuration of a single isolated PE-macroion complex. 
Two such complexes are then linked one by one to form a 
fiber with the help of straight linker chains and the only 
degree of freedom in this case is the relative rotation of 
unit cells around the linker chain. In both cases, a variety 
of helical structures are obtained for the fiber including, in 
particular, the zig-zag patterns. There are however quali- 
tative differences between the two models. 

In unconstrained minimization, the PE conformation 
is free to change on the sphere, i.e., upon changing the 
salt concentration or the macroion charge, the PE chain in 
each core complex may become increasingly more wrapped 
or dewrapped. As a result, the entry-exit angle of the PE 
strand varies accordingly, leading to dramatic changes in 
the overall fiber structure. As shown, the optimal fiber 
structure shows beads-on-a-string pattern or thick and 
loose solenoidal patterns at low salt. At about the in- 
verse screening length k = 0.6 nm -1 , the chain is already 
wrapped around the sphere for almost a complete turn. 
In this case, the spheres are aligned and form a fiber of 
diameter about 10 nm, which to some extent resembles 
the swollen 10 nm beads-on-a-string chromatin fiber. For 
salt concentrations within the physiological regime, zig- 
zag patterns are found as optimal structures and the fiber 



diameter is about 30 nm; this diameter and the fact that 
the fiber structure is zig-zag in the physiological regime 
thus appear to agree with those proposed for chromatin 
in Refs. [2112511251I271I251I25] . although the present model 
predicts a qualitatively different core particle structure. 
This is because the linker histone (absent in the present 
model) plays a key role in the compaction of chromatin 
into a 30 nm fiber by gluing together the entering and 
exiting strands of the DNA and thus sealing off the inter- 
nal structure of the DNA in the nucleosome core parti- 
cle. In our unconstrained model, the entering and exiting 
strands of the chain are displaced and divergent and thus, 
the density of spheres along the fiber appears to be much 
smaller (i.e., about one sphere per 11 nm along the fiber 
axis) than the value proposed in the cross-linker models 
(about 5-6 spheres per 11 nm [1511211261 12711251125 ] ) . For 
higher salt concentrations, the chain wrapping degree still 
increases and the spheres become highly packed along the 
fiber backbone (for fixed chain length per sphere) , leading 
to compact solenoidal structures. 

We have also studied the role of sphere charge and 
chain length per sphere and their influence on the fiber 
structure. An important result is that when the chain 
(linker) length per macroion is treated as a free param- 
eter, the binding energy (per sphere or unit cell) is found 
to take a non-convex shape at intermediate salt concen- 
trations when plotted as a function of the chain length 
per sphere. A simple common-tangent construction indi- 
cates a gas-liquid-type "phase coexistence" along the fiber, 
i.e., a part of the PE chain forms a dense complex with 
macroions phase-separating from a dilute phase along the 
fiber. 

In the constrained optimization model, the core parti- 
cle structure is fixed and only the relative rotation of the 
core particles (unit cells) is allowed in order to minimize 
the fiber Hamiltonian. As shown, for large sphere charge 
and within the physiological salt regime, the fiber shows 
zig-zag pattern and sphere densities of about 5-6 spheres 
per 11 nm length along the fiber axis. Indeed, both the 
global structure as well as the fiber density and diameter 
in this case resemble more closely the structures found or 
proposed for chromatin in Refs. [2112611271125 1125] . Inter- 
estingly, the dense zig-zag pattern results from the fact 
that the entering and exiting strands form a cross pattern 
leading to a small entry-exit angle; the form of the entry- 
exit strands in this case tends to resemble the stem-like 
form expected in chromatin in the presence of the linker 
histone [27, 28] and is therefore found to be a similarly im- 
portant factor in order to obtain dense fiber structures in 
our model. The constrained model is thus similar to the 
two-angle model proposed by Woodcock et al. [27J, ex- 
cept that here we employ a realistic interaction potential 
between the core particles. 

We emphasize again that the present results are ob- 
tained for a generic charged chain-sphere complex fiber 
and no direct comparison with the 30 nm chromatin fiber 
(which involves additional structural and specific details) 
should be attempted. Especially, the linker histone effects 
or other factors such as histone tails or the specific shape 
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of the histone octamer are not incorporated within the 
present model, and the constrained optimization model 
is only a simple way to mimic the linker-histone stabi- 
lization of the nucleosome. On the other hand, our model 
and its predictions are applicable to the mixtures of DNA 
and synthetic oppositely charged spheres and suggest that 
chromatin-like structures should also be observable in such 
systems. The main advantage of the present approach is 
that it accounts for the detailed wrapping state of the PE 
chain within a fiber as well as the essential intra-fiber in- 
teractions, and that it can easily be generalized to more 
sophisticated models. 

In this work, we have also neglected the excluded- 
volume interaction between chain segments. (Note that 
the chain-sphere volume interaction and the finite radius 
of DNA is taken into account via the hard-core radius 
of the sphere as explained in the text.) The chain-chain 
excluded-volume interaction is expected to play a role 
only at high salt concentration (typically beyond lOOmM 
monovalent salt for intermediate sphere charge) , where the 
chain wraps around the spheres more than two turns. (For 
a complete analysis and possible improvements of the ba- 
sic chain-sphere scheme used here see the discussions in 
Refs. [5l[5Sl[56].) 

Another interesting problem is to study the response 
of the fiber to an externally applied stress. Recent exper- 
iments have shown that the response of chromatin to an 
external mechanical stress can provide more insight into 
the DNA-histone interactions and the detailed conforma- 
tion of the fiber [1511^12211551151] . This problem has been 
considered in several recent theoretical studies [TTlll"B,36, 
38,54 . Finally, one may also consider an inhomogeneous 
fiber within the present model by assuming that the linker 
length may vary from one unit cell to the other. In this 
case, the fiber may exhibit overall deformations and large- 
scale conformational fluctuations [27] , 
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